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Abstract 

In  probabilistic  design  of  materials  for  fracture-critical  components  in  modem  military  turbine 
engines,  a  typical  design  target  risk  (DTR)  is  5x1 0~8  component  failures/engine  flight  hour.  This 
goal  highlights  the  essential  role  of  safety  in  a  design  process  that  simultaneously  strives  for  a 
balance  of  performance,  efficiency,  safety,  reliability,  and  affordability  over  the  entire  life  cycle 
of  the  engine.  Historically,  to  meet  safety  and  reliability  requirements,  the  design  and  life 
management  approaches  for  engine  materials  have  involved  large  experimental  programs  from 
which  statistically  based  life  limits  are  derived  by  extrapolation  from  mean  fatigue  behavior. 
However,  we  have  recently  found  that  the  fatigue  lifetime  under  a  given  test  condition  often 
exhibits  a  bimodal  form,  and  the  trends  in  mean  vs.  minimum  fatigue  lifetime  typically  respond 
differently  to  loading  and  microstructural  variables.  Under  such  circumstances,  the  underlying 
life-limiting  mechanisms  appear  to  follow  a  probabilistic  microstructural  hierarchy  that  is 
controlled  by  susceptibility  to  early  damage  and  growth  of  small  cracks.  These  findings  suggest 
that  significant  opportunities  may  exist  for  substantial  reductions  in  uncertainty  in  materials  life 
prediction  and  life  management.  The  potential  implications  of  the  new  findings  are  explored, 
and  a  number  of  possible  approaches  are  suggested  for  incorporating  the  insights  of  life-limiting 
fatigue  into  methods  of  integrated  computational  materials  engineering  (ICME)  that  support 
optimized  life-cycle  design  of  materials  and  components  in  turbine  engines.  Potential  benefits  of 
this  approach  appear  to  include  substantial  improvements  in  model  accuracy,  coupled  with 
reduced  requirements  for  materials  testing,  potentially  leading  to  a  significant  reduction  in  the 
time  and  cost  to  develop,  validate,  transition,  and  implement  new,  more  fatigue,  resistant  alloys. 
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Introduction 

A  bright  future  for  materials  design 

Looking  to  the  future  of  turbine  engine  materials,  it  appears  that  there  are  many  potential 
breakthroughs  in  the  underlying  science  and  technology  of  materials  design  that  could  offer 
major  opportunities  for  improvement  beyond  today’s  engineering  state-of-the-art.  For  example, 
future  designs  are  likely  to  benefit  substantially  from  a  wave  of  anticipated  advancements  in 
probabilistic  modeling  and  simulation,  micro-experimentation,  high  fidelity  materials 
characterization,  mechanism-based  modeling,  and  information  science  and  technology,  all  of 
which  will  be  spurred  on  by  a  number  of  major  new  initiatives  in  science  and  technology[l-4]. 
Through  unprecedented  levels  of  model  integration  between  the  science  and  engineering 
disciplines,  it  should  be  possible  to  design  future  propulsion  systems  to  be  more  adaptable  to  new 
usage  requirements,  to  have  higher  fuel  efficiency,  to  be  inherently  more  reliable  and  verifiable, 
to  be  manufactured  and  certified  more  quickly,  and  to  require  far  less  logistics  support  over  the 
life  cycle. 
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In  addition,  real-time  tracking  of  engine  performance,  flight  profdes,  and  diagnostics  of 
structural/materials  health  are  expected  to  give  operators  and  logisticians  unsurpassed  insights 
for  management  and  sustainment  of  individual  assets,  as  part  of  a  comprehensive  philosophy  of 
condition  based  management  +  prognosis  [5-9],  including  the  concept  of  a  digital  twin  to 
represent  individual  assets  [10].  Undergirding  these  advancements,  the  burgeoning  field  of 
Integrated  Computational  Materials  Engineering  (ICME)  [11,  12]  will  enable  game-changing 
materials-design  concepts,  including  the  ability  to  tailor  location-specific  material  properties  that 
are  optimized  to  meet  component-specific  life-cycle  design  requirements.  Perhaps  the  greatest 
underlying  challenge  for  the  science  and  engineering  community,  however,  will  be  to  achieve  the 
full  potential  of  a  material’s  capabilities  along  with  simultaneous  requirements  for  safety, 
reliability,  and  affordability  over  the  full  life  cycle.  In  fact,  the  growing  importance  of 
affordability  was  recently  underscored  in  the  U.S.  Air  Force  forecast  Technology  Horizons, 
Essential  Focus  Areas  for  Air  Force  S&T  Investment :  “. .  .Air  Force  science  and  technology  over 
the  next  decade  will  need  to  focus  as  much  on  advancing  technologies  that  enable  reduced  Air 
Force  operating  costs  as  on  technologies  that  support  more  traditional  development  of  new 
systems  or  capabilities."  [13] 

To  help  meet  this  challenge,  the  engine  structural  design  community  is  immersed  in 
building,  improving,  and  exploiting  comprehensive  probabilistic  design  tools  for  engine 
components  [e.g.,  [14]].  Underpinning  such  computer  codes  are  the  variabilities  and 
uncertainties  in  the  input  values,  including  probabilistic  descriptions  of  engine  missions/loads, 
aerodynamic  flow,  thermal  and  stress  analysis,  combustion  effects,  and  materials  performance  in 
harsh  environments.  At  the  structures  and  system  level,  new  concepts  for  robust  design  [15]  and 
prognosis  [5]  are  emerging.  In  this  context,  Bayesian  updating  offers  a  rigorous  approach  to 
refine  model  accuracy  based  on  continuous  feedback  from  diagnostic,  maintenance,  and 
operational  data  [16].  Asa  key  part  of  these  improvements,  better  models  of  materials 
performance  and  uncertainty  are  also  needed. 

Case  study:  turbine  engine  materials  -  building  on  a  history  of  engineering  design  and  operation 

Life-cycle  management  procedures  used  for  today’s  advanced  turbine  engines  have 
achieved  a  superb  record  of  safety  and  reliability,  and  the  underlying  models  and  methods  have 
been  continually  improved  over  decades  of  design  and  operation.  This  process  has  culminated  in 
highly  validated  engineering  standards  of  practice,  supported  by  large  databases  of  material 
properties.  Consequently,  these  engineering  practices  are  highly  tolerant  to  variations  and 
uncertainties  in  materials  and  structures,  and  catastrophic  failures  of  components  have  been 
virtually  eliminated.  At  the  same  time,  owing  largely  to  the  mature  age  of  some  of  the  fleets, 
engine-system  life-cycle  costs  are  growing  at  an  annual  rate  that  is  well  above  inflation,  and  this 
cost  burden  threatens  to  stifle  efforts  to  recapitalize  and  modernize  future  engine  fleets  [17]. 
Given  such  affordability  challenges,  there  is  a  growing  need  to  reduce  life-cycle  costs.  However, 
history  indicates  that  as  much  as  80%  of  an  engine’s  life  cycle  cost  is  established  during  the 
design  and  early  development  of  the  engine,  and  the  cost  to  correct  design  defects  afterward 
becomes  exponentially  greater  as  development  proceeds  [17]. 

Within  the  U.S.  Department  of  Defense,  the  decades  of  engine  design  and  engineering 
experience  have  molded  the  specifications  known  as  the  Propulsion  System  Integrity  Program 
(PSIP)  [18]  and  its  associated  document  [19],  which  guide  the  life-cycle  design  and  management 
of  turbine  engines.  The  PSIP  design,  depicted  schematically  in  Figure  1,  is  the  beginning  of  a 
comprehensive  life-cycle  design  and  management  process  that  includes:  (i)  Design  Information 
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and  Development  Planning;  (ii)  Design  Characterization  and  Analysis;  (iii)  Components  and 
Subsystems  Testing;  (iv)  Propulsion  System  Ground  and  Flight  Testing,  and  (v)  the  Engine  Life 
Management  Plan.  Engine-material  development,  selection,  characterization,  modeling,  and 
validation  are  concentrated  in  Tasks  (i)  and  (ii),  and  first-tier  material  properties,  such  as 
strength,  stiffness,  and  density,  tend  to  drive  initial  design.  Second-tier  material  properties,  such 
as  cycle-,  time-,  temperature-,  and  environmental-dependence,  are  prominent  factors  throughout 
the  design  and  life  cycle,  and  these  properties  typically  define  engine  overhaul  intervals,  affect 
the  average  time  on  wing,  and  dictate  the  ultimate  lifetimes  of  components. 

For  safety-  or  mission-critical  components,  such  as  rotors  and  rotor  spacers,  the  PSIP 
design  life  is  driven  largely  by  major  load  cycles  produced  by  variations  in  engine  rotational 
speed,  wherein  a  ground-air-ground  flight  is  counted  as  1  cycle,  and  minor  throttle  excursions 
during  flight  count  a  partial  cycles.  Such  cycles  may  be  summarized  as  total  accumulated  cycles 
(TACs),  and  a  typical  flight  produces  approximately  2-3  TACs/hour.  To  satisfy  performance 
and  safety  requirements,  a  design  life  for  major  rotating  components  in  high  performance 
engines  is  nominally  8000  TACs,  with  a  routine  mechanical  overhaul  of  the  engine  prescribed  at 
4000  TACs.  Given  these  target  lifetimes,  the  associated  maximum  stresses  in  rotor  components 
tend  to  approach  the  material  yield  strength  at  key  locations,  such  as  boltholes  and  other  stress 
concentrators,  and  relatively  high  maximum  stresses  occur  throughout  the  components. 

With  regard  to  the  materials  life  cycle,  it  is  well  known  that  advanced  Ti  and  Ni  alloys  in 
turbine  engine  rotors  exhibit  a  substantial  variation  in  their  long-term  cyclic-  and  time-dependent 
properties,  such  as  resistance  to  fatigue  damage  and,  to  a  lesser  extent,  resistance  to  creep 
deformation.  The  PSIP  specifications  deal  with  such  variation  by  managing  to  a  criterion  based 
on  statistical  lower-bound  behavior.  For  example,  for  design  of  safety-  or  mission-critical 
components,  these  specifications  impose  simultaneous  requirements  for  (i)  a  statistical  minimum 
safe-life  fatigue  and  (ii)  damage  tolerance  to  small  surface  or  comer  cracks.  Under  the  safe-life 
criterion,  the  minimum  LCF  life  is  defined  as  -3  standard  deviations  from  the  mean  value  of  low 
cycle  fatigue  lifetime,  or  alternatively,  the  1/1000  (B0.1)  lower-bound  capability,  based  on  a 
database  of  tests  over  the  range  of  engine  design  conditions.  Thus,  the  PSIP  specification 
virtually  guarantees  material  safety  and  reliability  under  low  cycle  fatigue,  although  the  reliance 
on  large  statistical  databases  can  be  costly  and  time  consuming  to  produce. 

Engines  are  designed  to  meet  a  target  mission  profile  of  loads,  but  actual  mission  usage 
invariably  change  as  the  operational  needs  change,  making  tracking  of  actual  mission  usage  an 
essential  input  to  the  life  analysis.  Future  engines  are  expected  to  accommodate  this  need  with 
sophisticated  usage  tracking,  and  the  external  forces  and  boundary  conditions  on  engine 
components  are  expected  to  be  well  defined,  given  the  highly  evolved  state  of  mechanical  design 
and  stress  analysis.  However,  the  understanding  and  control  of  variability  in  materials  response 
to  such  loading  is  much  less  mature. 

In  the  quest  for  improved  models  to  address  the  goals  of  safety,  reliability,  and 
affordability,  safety  is  likely  to  remain  as  the  number  one  requirement  for  manned  aircraft.  For 
example,  as  noted  in  the  PSIP  specification,  the  limit  for  the  design  target  risk  (DTR)  for 
probabilistic  design  of  a  component  “was  set  at  a  value  of  5x1  O'8  predicted  occurrences  per  EFH 
[Engine  Flight  Hour]  for  a  recent,  advanced  single-engine  fighter.”  Similarly,  specifications  for 
commercial  engine  components  set  the  minimum  DTR  at  lxl 0'9/flight  [20].  Such  rigorous 
safety  requirements  underscore  the  importance  of  fatigue  life  limits  in  design  and  operation  of 
turbine  engines. 
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Life  limits  and  the  duality  of  fatigue 

Under  previous  efforts  on  life-limiting  material  behavior,  we  have  found  that  the 
variability  in  fatigue  response  of  engine  alloys  at  a  given  stress  level  often  tends  to  separate  into 
bimodal  distributions  representing  minimum  behavior  and  a  mean-life  dominated  behavior, 
termed  Type  I  and  Type  II  behavior,  respectively.  Thus,  the  probably  density  representing  a 
material’s  safe-life  fatigue  capability  shown  in  Fig.  1  was  often  composed  of  two  underlying 
probability  densities,  as  depicted  schematically  in  Fig.  2.  This  duality  in  fatigue,  also  known  as 
bimodal  or  competing-modes  fatigue,  has  been  observed  in  a  wide  range  of  alloys,  including  the 
superalloys:  Rene’95  [21,  22],  Rene’88DT  [23],  IN100  [8,  24-29],  Waspaloy  [30],  the  single 
crystal  alloy  PWA  1484  [31],  the  titanium  alloys  Ti-10-2-3  [23,  32-34],  Ti-6Al-2Sn-4Zr-6Mo 
[35-42],  Ti-6A1-4V  [43-47]  gamma  titanium  aluminides  [48,  49],  and  the  aluminum  alloy  7075- 
T651  [45],  and  others.  Although  such  a  separation  of  fatigue  response  has  been  known  for  some 
time,  the  significance  of  this  behavior  has  not  yet  been  generally  captured  in  the  strategies  for 
fatigue  design  of  turbine  engine  materials. 

In  exploring  this  behavior  in  turbine  engine  alloys,  we  have  developed  a  physically-based 
approach  for  describing  fatigue  variability,  and  this  approach  has  been  integrated  into  a  Monte 
Carlo  probabilistic  life  prediction  model  for  materials[35,  40,  50].  The  primary  tenets  of  our 
fatigue  variability  description  are  summarized  by[40]:  (i)  under  nominal  microstructural  and 
loading  conditions,  a  hierarchy  of  local  deformation  heterogeneities  develop  in  a  fatigue  sample, 
corresponding  to  certain  microstructural  arrangements  or  features;  (ii)  as  a  result  of  this 
hierarchy,  a  probability  exists  of  an  extreme  microstructural  arrangement  that  may  initiate  a 
crack-growth  dominated  failure,  producing  a  life-limiting  failure  distribution  for  the  material; 
and  (iii)  there  are  often  different  degrees  of  influence  of  a  given  variable  on  minimum  and  the 
mean  lifetimes.  For  example,  the  effect  of  stress  level  on  mean  vs.  life-limiting  fatigue  typically 
produces  quite  different  responses,  and  such  differing  response  of  minimum  and  mean  fatigue 
behavior  has  been  observed  for  other  important  variables,  including  microstructure,  surface 
treatment,  temperature,  and  dwell  time.  This  understanding,  when  applied  to  fatigue  design  of 
materials,  suggests  that  compositional,  processing,  surface  treatment,  and  microstructural 
modifications  should  consider  the  lower-lifetime  limit,  and  that  the  life-limiting  behavior  may 
not  be  accurately  represented  by  purely  statistical,  mean-fatigue-based,  approaches. 

The  goal  of  the  current  paper  is  to  (i)  overview  general  findings  on  fatigue  life  limits  of 
alloys  used  in  turbine-engine  rotor  components,  (ii)  to  highlight  the  effects  and  implications  of 
such  fatigue  life  limits  for  a  specific  high  strength  titanium  alloy,  and  (iii)  to  identify  possible 
avenues  for  improvement  in  materials-model  accuracy.  Thus,  the  focus  will  be  on  challenges 
and  opportunities  to  reduce  variability  in  performance,  and  uncertainty  of  predictions,  of  fatigue 
life-cycle  capability  of  turbine-engine  rotor  materials,  consistent  with  the  demands  of  a  high- 
performance  single-engine  aircraft. 

Material  and  experimental  methods 

The  focus  of  this  investigation  was  a  pancake  forging  of  alloy  Ti-6Al-2Sn-4Zr-6Mo  processed 
and  heat-treated  to  a  duplex  a+p  phase  microstructural  condition,  as  shown  in  Fig.  3.  The 
micro  structure  consisted  of  primary,  equiaxed  a  particles  (ap),  and  lamellar  colonies  of  fine  a/p 
laths,  which  transformed  from  prior-p  grains  during  heat  treatment.  The  average  size  of  the  ap 
was  about  4  pm  and  that  of  the  colonies  was  about  15  pm.  The  volume  fractions  of  ap  and  the 
colony  phases  were  approximately  0.3  and  0.7,  respectively.  This  produced  a  bulk  yield  strength 
of  1040  MPa  and  an  ultimate  strength  of  1250  MPa. 
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Closely  spaced,  cylindrical  specimen  blanks  of  this  material  were  machined  from  the 
forging,  and  sections  of  Ti-6A1-4V  rod  were  inertia  welded  to  the  ends  as  extensions  for 
gripping.  Subsequently,  the  inertia-welded  blanks  were  machined  and  finished  by  low-stress 
grinding  into  either  cylindrical  dog-bone  specimens  for  fatigue  testing,  or  similar  small-crack 
specimens,  which  were  fabricated  with  surface  flats  to  accommodate  crack  growth  measurements 
via  surface  replication  using  standard  acetate  replication  procedures.  Surface  residual  stresses 
from  machining  were  removed  by  electropolishing  the  machined  specimens  to  remove  a  50  pm 
thick  surface  layer  on  the  specimen  gage  sections.  The  small-crack  tests  were  performed  using 
either  natural  initiation  or  by  initiation  from  multiple  small  surface  notches,  which  had  been 
machined  into  the  surface  using  a  focused  ion  beam  (FIB)  to  create  residual-stress-free  crack 
starters. 

To  assess  the  inherent  effect  of  the  Ti-6246  alloy  microstructure  on  the  variability  in 
fatigue  response,  a  series  of  tension-tension  fatigue  tests  was  performed  at  each  of  nine  stress 
levels  and  a  stress  ratio  of  min/maximum  stress  of  0. 1 .  The  specific  maximum  stresses  for  these 
tests  were  1040,  925,  900,  860,  820,  700,  650,  600,  and  550  MPa.  The  test  frequency  was  20  Hz 
for  maximum  stresses  above  750  MPa,  while  tests  below  this  stress  were  performed  in  ultrasonic 
resonance  at  20  KHz.  At  the  highest  stresses,  surface-connected  a  particles  typically  served  as 
the  primary  sites  for  crack  initiation  (e.g.,  Fig.  4),  while  some  of  the  specimens  failed  from 
internal  initiation  sites  under  maximum  stresses  of  860  MPa  or  lower.  Further  details  on  the 
microstructure,  experimental  methods,  results,  and  analysis  are  provided  in  [37,  39-41,  51], 
which  served  as  the  backbone  of  the  results  to  be  presented  below. 

Results 

Figure  5  presents  an  S-N  plot  of  the  resulting  fatigue  lifetimes  for  all  stress  levels  (amax  of 
550  to  1040  MPa)  and  both  frequencies.  As  shown,  the  variability  in  fatigue  lifetime  exceeded 
three  orders  of  magnitude  at  the  intermediate  stress  levels,  while  the  variability  in  lifetime  at  the 
highest  stress  level  (1040  MPa)  was  extremely  limited.  The  increase  in  variability  with  reducing 
stress  has  been  observed  widely  in  the  literature,  although  the  extent  of  the  increase  is  likely 
maximized  for  these  tests,  due  to  the  residual-stress-free  electropolished  surfaces.  This  surface 
condition  was  chosen,  however,  precisely  because  it  eliminated  effects  of  surface  irregularities 
and  residual  stresses  and  helped  to  reveal  the  inherent  effects  of  alloy  microstructure.  The  nature 
of  the  statistical  variations  in  the  fatigue  lifetimes  at  each  stress  level  is  elucidated  in  the 
probability  plot  of  Fig.  6.  As  shown,  the  lifetimes  from  the  highest  stresses  exhibited  a  straight- 
line  trend  on  the  lognormal  probability  scale,  while  bimodal  trends  of  the  data  tended  to  emerge 
at  intermediate  stress  levels.  As  shown  on  an  expanded  scale  in  Fig.  7,  the  bimodal  behavior  was 
most  evident  for  the  fatigue  datasets  at  amax  =  860  and  820  MPa.  For  each  of  these  datasets, 
three  regions  of  fatigue  behavior  are  distinct  (Type  I,  transition,  and  Type  II),  with  the  ratio  of 
the  data  points  in  the  Type  I  vs.  Type  II  regions  being  stress-dependent. 

These  separate  regions  of  fatigue  response  were  use  as  a  basis  for  isolating  the  short-  and 
long-life  trends  alone,  as  re-plotted  in  Fig.  8.  As  shown,  when  the  Type  I  and  Type  II  regions 
are  plotted  separately,  the  data  appear  nominally  straight  on  the  lognormal  probability  scale, 
indicating  that  each  region  is  unimodal.  Note  that  the  data  from  the  amax  =  860  and  820  MPa 
conditions  are  plotted  twice,  since  both  Type  I  and  Type  II  regions  were  clearly  present  at  these 
stress  levels.  In  comparison  with  the  full  range  of  fatigue  variability  present  in  the  complete 
datasets  of  Fig.  6,  Fig.  8  shows  that  two  regions  of  behavior  tend  follow  distinctly  different 
trends,  appearing  to  represent  a  life-limiting  (Type  I)  behavior,  plus  a  final  period  of  fatigue 
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exhaustion  where  all  specimens  eventually  failed.  Interestingly,  in  terms  of  the  lognormal 
probability  scale  of  the  plot,  there  was  limited  variation  of  the  slopes  of  the  data  for  the  two 
regions  of  behavior  (life-limiting  and  fatigue-exhaustion),  and  the  slopes  of  the  two  regions  were 
remarkably  similar.  Although  all  specimens  failed  (no  run-out  specimens),  the  longest-lived 
specimens  may  have  been  approaching  a  fatigue  limit. 

Figure  9  shows  the  data  from  the  820  MPa  tests,  and  includes  a  lognormal  probability 
regression  line  with  95%  confidence  intervals  on  the  regression.  Clearly,  such  a  regression  is 
inappropriate  for  bimodal  data.  However,  materials  and  component  designers  rarely  if  ever  have 
the  luxury  of  having  so  many  repeated  tests  at  an  individual  fatigue  condition.  With  fewer  test 
results,  the  bimodal  behavior  would  probably  not  be  evident,  and  a  regression  similar  to  the  one 
shown  would  likely  appear  reasonable.  In  this  case,  one  can  see  that  extrapolation  to  the  B0.1 
probability  of  failure  would  give  approximately  3500  cycles,  while  the  corresponding  lower 
bound  of  the  95%  confidence  interval  was  about  721  cycles.  Alternatively,  taking  maximum 
benefit  of  the  large  number  of  data  points  at  this  condition  allows  one  to  focus  on  the  single, 
short-life  (Type  I)  distribution.  Figure  10  shows  a  probability  plot  of  the  6  shortest-life  data 
alone,  and  in  this  instance  the  estimated  B0.1  life  was  greater  than  104  cycles,  and  the  minimum 
95%  confidence  estimate  was  approximately  5660  cycles,  or  about  8  times  the  corresponding 
value  indicated  by  the  unimodal  analysis  of  the  full  data  set.  Thus,  this  example  implies  that 
existence  of  bimodal  fatigue  behavior  may  indicate  that  the  material’s  inherent  minimum  lifetime 
is  substantially  greater  than  otherwise  thought. 

Analysis 

Modeling  bimodal  fatigue 

A  key  opportunity  in  the  design  of  future  materials  is  to  understand,  isolate,  and  model  the 
sources  of  variability  in  fatigue  damage  and  lifetime.  As  noted  earlier,  we  previously  established 
a  general  computational  framework  for  physics-based  understanding  of  fatigue  variability  [24, 
40],  In  this  approach,  the  lifetime  distribution  was  modeled  as  a  superposition  of  the  probability 
densities  of  the  crack-growth  lifetime  and  the  mean  lifetime.  The  total  lifetime  density  was 
represented  by: 


ft(N)  =  Vlfl(N)  +  Pmfm(N)  (1) 

where,  f(N),  f(N),  and  fm(N)  are  probability  densities  representing  the  total  lifetime,  crack  growth 
lifetime,  and  the  mean  lifetime,  respectively.  The  factors,  pi  and  pm  are  the  probabilities  of 
occurrence  of  the  respective  responses.  Salient  features  of  the  model  are  (i)  that  the  effect  of  any 
variable  on  lifetime  distribution  is  separated  into  its  influence  on  small-plus-large  crack  growth 
vs.  the  mean  lifetime  behaviors,  and  (ii)  that  the  lifetime  distribution  is  modeled  as  a 
superposition  of  these  two  responses.  From  a  design  life  perspective  (for  example,  a  1/1000, 
B0.1,  lifetime),  it  can  be  shown  that  the  prediction  of  life  limit  by  the  crack  growth  lifetime 
density,  f(N)  (i.e.,  in  the  limit  of  pi  — »  1)  forms  a  lower  bound  of  predictions  by  the  bimodal 
model  [35],  as  discussed  below. 

It  is  well  known  that  the  total  fatigue  life  of  a  specimen  is  composed  of  the  sum  of  the  cycles 
from  the  successive  stages  of  damage  initiation  and  crack  growth  to  failure,  such  that  the  total 
fatigue  life  (Nt)  is  given  by: 


Nt  —  Ni  +  Nsc  +  Nlc 
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where  N„  Nsc,  and  Nlc  are  the  cycles  spent  in  crack  initiation,  small  crack  growth,  and  large- 
crack  growth,  respectively.  It  is  also  recognized  that  Nt  may  be  further  subdivided  into  a  crack- 
incubation  plus  a  crack-nucleation  phase  in  some  instances,  and  that  small  cracks  may  be  defined 
as  being  microstructurally,  physically,  mechanically,  and/or  chemically  small  [52-54],  For  a 
given  material  and  loading  condition,  the  variability  in  total  fatigue  lifetime,  Vt,  is  a  sum  of  the 
variabilities  of  the  individual  stages  of  damage  and  cracking,  and  variability  of  the  individual 
terms  in  Eq.  1  tends  to  reduce  from  left  to  right,  such  that  V(N,)  >  V(Nsc)  >  V(Nlc)- 

As  noted  earlier,  for  a  typical  8000-cycle  design  life  of  turbine-engine  rotor  materials,  the 
corresponding  maximum  stresses  approach  the  material’s  yield  stress.  Under  these 
circumstances,  it  appears  that  the  number  of  cycles  to  crack  initiation  from  a  microstructural- 
scale  feature  may  occasionally  approach  zero  for  the  life-limiting  fatigue.  Assuming  that  A  ~  0, 
one  can  calculate  the  expected  lifetime  by  integrating  the  equation 

,,  A  J  da  7  da 

Nt  —  0  +  I - h  I  -  /'j'v 

where  the  two  integrals  represent  the  life  from  small-crack  and  large-crack  growth,  respectively, 
at,  is  the  crack  initiation  depth  from  a  microstructural  or  mechanical  defect,  asc  is  the  transition 
crack  size  between  small-  and  large-crack  growth,  and  fsc  and  fie  are  the  crack  growth  functions 
for  the  small-  and  large-crack  size  ranges,  respectively.  For  the  Monte  Carlo  model  of  the  life- 
limiting  fatigue  process,  the  distribution  of  crack  initiation  sizes  (a,)  was  represented  by  the  red 
curve  of  Fig.  1 1.  A  statistical  representation  of  the  small-crack  growth  behavior  shown  in  Fig. 

12  was  obtained  using  maximum  likelihood  methods  to  obtain  a  correlated  Paris-law  function  to 
represent  the  variability  in  crack  growth  behavior,  while  the  large-crack  behavior  was  treated  as 
deterministic  (mean  value).  For  each  stress  level,  at  least  10,000  Monte  Carlo  simulations  were 
used  to  calculate  a  distribution  of  expected  minimum  total  lifetimes,  assuming  N,  =  0.  Figure  13 
plots  these  predicted  minimum  lifetimes  in  comparison  with  the  omax  =  820  MPa  life-limit  data 
of  Fig.  10.  As  shown,  the  probabilistic  prediction  agrees  extremely  well  with  the  life-limit  data, 
which  tends  to  corroborate  the  assumption  that  A,  was  at  least  approaching  zero  for  these  tests. 

Figure  14  re-plots  the  S-N  data  of  Fig.  5,  and  shows  a  schematic  of  bimodal  regressions 
of  the  fatigue  behavior  for  each  stress  level.  The  three  lines  at  the  left  of  the  plot  represent 
Monte  Carlo  model  predictions  of  the  0.1%,  10%,  and  50%  probabilities  of  failure,  respectively, 
from  left  to  right.  As  shown,  all  of  the  fatigue  data  are  bounded  approximately  by  the  P(failure) 

=  10%  line.  In  addition,  the  slopes  of  the  Monte  Carlo  predictions  vs.  stress  level  agree  well 
with  the  slope  of  the  Type  I  fatigue  lifetime  line  representing  the  mean  behavior  of  the  life- 
limiting  populations  alone.  Among  the  implications  of  the  steeper  slope  for  life-limiting  fatigue 
is  a  reduced  sensitivity  to  stress  level,  as  compared  to  a  traditional  view  of  the  mean  S-N  fatigue 
behavior. 

Also  as  shown  in  the  figure,  the  proportions  of  data  falling  in  the  short-  vs.  long-life 
regions  of  S-N  behavior  shift  toward  longer  lifetimes  as  stress  is  reduced.  Recalling  Eq.  1 ,  this  is 
equivalent  to  reducing  the  probability  (pi)  of  a  given  fatigue  data  point  being  in  the  short  life 
region.  Based  on  a  nonlinear  bimodal  regression  of  the  fatigue  data  at  each  stress  level,  an 
estimate  of  the  dependence  of  pi  on  stress  level  is  plotted  in  Fig.  15.  As  shown,  the  maximum 
stress  at  which  the  short-  and  long-life  behaviors  become  equally  likely  is  approximately  820 
MPa,  which  is  about  79%  of  the  yield  strength  for  the  Ti-6246  alloy  tested. 
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Figure  16a  presents  the  results  of  1000  Monte  Carlo  simulations  of  crack  growth 
behavior  at  820  MPa,  and  Fig.  16b  shows  an  expanded  region  of  the  same  plot.  As  shown,  for 
the  distribution  of  crack  initiation  depths  (up  to  approximately  15  pm  in  this  case),  the  minimum 
lifetime  is  more  than  14000  cycles,  which  is  well  in  excess  of  the  typical  8000  cycle  (TAC) 
design  lifetime.  Recognizing  that  the  estimate  of  pi  at  this  stress  level  is  0.5,  half  of  all 
specimens  tested  at  this  condition  would  be  expected  to  lie  within  the  short-life  distribution. 

Since  traditional  low-cycle-fatigue  safe-life  design  targets  the  1/1000  statistical  minimum  for 
initiation  of  a  0.89  mm  deep  crack,  this  should  correspond  to  the  minimum  values  of  the  Monte 
Carlo  predictions.  As  such,  an  alternative  interpretation  of  the  1/1000,  or  -3  standard  deviation, 
LCF  safe-life  criterion  is  that  for  these  cases,  in  all  likelihood,  beginning  at  the  size  of  a 
microstructural  feature  (~  5-15  pm),  a  crack  would  have  been  growing  for  most  of  the  safe-life 
period.  There  is  a  notable  caveat  in  this  regard,  however,  since  the  current  results  are  for  a 
stress-free  surface  condition,  while  actual  components  typically  have  highly  compressive  surface 
residual  stresses  from  machining  plus  shot  peening,  as  discussed  below. 

Effects  of  surface  residual  stresses  produced  by  shot  peening 

Surface  treatments  such  as  shot  peening  are  used  widely  to  retard  crack  initiation  and 
propagation  in  the  surface  region,  and  thereby  increase  the  lifetime.  However,  our  findings  from 
other  research  suggest  that  the  effect  of  shot  peening  on  the  lifetime  distribution  should  be 
represented  in  terms  of  its  effect  on  the  mean  and  the  life-limiting  behavior  [42]  .  In  work  on  the 
same  Ti-6246  as  reported  in  the  current  paper,  shot-peening  increased  both  the  minimum  and 
mean  lifetimes  by  suppressing  crack-initiation  in  the  surface  region  and  forcing  crack-initiation 
to  the  subsurface  in  many  samples  [42],  Although  failure  by  the  subsurface  mechanism 
improved  the  mean  response,  the  life-limiting  failures  under  shot-peening  occurred  by  a  different 
surface-crack-initiation  mechanism  in  which  cracking  was  found  to  originate  as  an  interaction 
between  shot-peening  and  the  local  microstructure,  which  created  local  regions  of  surface  micro¬ 
cracking.  These  surface  crack  initiation  sizes  appeared  to  be  significantly  larger  than  the  average 
primary  a  grain  size  (~5  pm)  that  controlled  crack  initiation  in  the  residual  stress  free  condition. 
The  limiting  lifetimes  under  shot  peening,  therefore,  approached  to  within  approximately  1.7 
times  longer  than  those  under  the  residual-stress  free  condition.  In  other  words,  a  portion  of 
beneficial  effect  of  surface  residual  stresses  from  shot  peening  was  partially  offset  by  the 
creation  of  surface  cracks  from  shot  peening  damage.  The  primary  benefit  of  the  residual  stress 
appeared  to  be  due  to  the  clamping  stresses  acting  perpendicular  to  the  crack  plane,  which 
reduced  the  effective  range  of  crack  opening  and  stress  intensity  factor  (AKeff). 

Opportunities  for  improvement:  sources  of  variability  and  uncertainty  in  material  life  limits 
Given  that  life  limits  play  such  a  crucial  role  in  the  design  and  management  of  turbine  engine 
rotor  materials,  it  appears  that  the  results  outlined  in  this  paper  represent  a  mere  first  step  in  the 
potential  to  reduce  and  model  the  variability  in  materials  fatigue  performance.  By  focusing  on 
the  mechanics  and  physics  that  control  the  short-life  tail  of  a  material’s  fatigue  capability,  it 
appears  possible  to  characterize  individual  contributions  to  life  limits  with  a  high  level  of 
fidelity,  potentially  leading  to  increases  in  capabilities  to  control  and  fully  utilize  inherent 
materials  properties.  Two  primary  avenues  for  further  development,  maturation,  and  transition 
of  the  concepts  of  life-limiting  fatigue  are:  (i)  through  the  integrated  design  and  certification  of 
materials  and  structures  for  future  systems  and  (ii)  through  improvements  to  methods  of  life- 
cycle  management  and  prognosis. 
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As  discussed  above,  beginning  with  the  residual-stress-free  case,  one  can  start  to  isolate 
effects  of  microstructure  on  inherent  fatigue  variability.  This  special  case  presents  the 
opportunity  to  identify  life-limiting  microstructural  neighborhoods  and  effects  of  microstructural 
gradients  on  fatigue.  As  shown,  the  inherent  life  limits  under  low  cycle  fatigue  appear  to  be 
governed  by  the  growth  of  a  dominant  crack  from  an  initial  crack  size  given  by  the 
microstructural  feature  that  controls  crack  formation.  Thus,  the  approach  of  segmenting  total 
lifetime  into  the  components  of  initiation,  small-,  and  large-crack  growth  offers  a  very  systematic 
process  for  characterizing  and  modeling  the  many  additional  effects  that  are  often  difficult  to 
quantify  independently.  Such  models  should  provide  new  interpretations  of  fatigue  behavior  that 
have  traditionally  been  based  primarily  on  extrapolations  from  mean- fatigue  properties.  This 
could  be  significant  in  some  instances,  since  minimum  and  mean  fatigue  responses  differ  based 
on  microstructural  and  loading  variables. 

In  light  of  the  rapidly  growing  science  and  technology  offered  by  Integrated  Computational 
Materials  Engineering,  the  modeling  of  life-limiting  fatigue  of  materials  presents  many 
opportunities  for  the  development,  transition,  and  implementation  of  new  understanding  for  both 
military  and  commercial  applications,  including: 

•  Quantitative  2-  and  3-dimensional  characterization  and  representation  of  life-limiting 
microstructural  features  vs.  nominal  and  off-nominal  microstructures 

•  Nondestructive  characterization  of  materials  deformation,  damage  states,  microstructural 
features,  and  residual  stresses 

•  High  resolution  microscale  testing  of  microstructural  features,  including  digital  image 
correlation  for  precision  characterization  of  crystallographic  mechanisms  of  deformation  and 
damage 

•  Crystal  plasticity  modeling  of  microstructural  features,  including  life-limiting  local 
neighborhoods  and  gradients 

•  Analysis  and  modeling  of  effects  of  surface  and  bulk  residual  stresses 

•  Effects  of  geometric  stress  concentrations,  such  as  bolt  holes  and  notches 

•  3 -dimensional  fracture  mechanics  of  complex  geometric  features 

•  Effects  of  complex  mission  loading  and  histories,  including  thermomechanical  fatigue  and 
dwell  hold  times 

•  Advanced  validation  materials  science,  including  selective  analysis  of  the  damage  states  of 
materials  from  retired  components  to  provide  model  feedback 

•  Information  integration  and  data  mining 

•  Methods  of  material  state  awareness  supporting  Condition  Based  Management  +  Prognosis 

•  Probabilistic  modeling  of  multi-scale  domains:  microstructure,  geometric  feature, 
component,  and  system 

•  Model-based  probabilistic  simulation  and  validation  of  materials  and  components  designs 

Conclusions 

This  paper  used  electropolished  specimens  of  the  high-strength  titanium  alloy  Ti-6A1- 
2Sn-4Zr-6Mo  to  explore  fundamental  variabilities  and  uncertainties  in  microstructurally  based 
fatigue  life  limits  under  stresses  and  lifetimes  representative  of  engine  rotor  design  (nominally 
8000  major  cycles,  (TACs)).  In  this  life  range,  it  is  well  known  that  the  total  fatigue  lifetime 
(Nt)  may  be  partitioned  as  (Eq.  2): 
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Nt  —  Ni  +  Nsc  +  Nlc 


where  Nh  Nsc,  and  Nlc  are,  respectively,  the  numbers  of  cycles  spent  in  crack  initiation,  small- 
crack  growth,  and  large-crack  growth  to  failure.  In  the  test  results  presented  here,  N  was  often 
the  dominant  term  contributing  to  the  mean  fatigue  lifetime,  but  the  fatigue  mean  behavior 
sometimes  also  encompassed  a  shorter-life  fatigue  response  wherein  Nt  was  small.  This 
minimum  vs.  mean  bimodal  fatigue  behavior  appeared  to  respond  differently  to  stress  level, 
which  may  be  expected  depending  on  the  relative  contribution  of  N  to  the  total  fatigue  lifetime. 
In  related  research,  we  have  also  observed  that  the  trends  in  minimum  vs.  mean  fatigue  response 
may  differ  as  a  function  of  microstructural  and  usage  variables. 

For  turbine-engine  rotor  materials,  traditional  design  practice  defines  the  safe-life  fatigue 
limit  as  a  statistical  minimum  life  to  form  a  small  surface  or  comer  crack.  Alternatively,  from  a 
materials-microstructure  perspective,  one  might  define  the  minimum  size  of  a  small  crack  as  the 
dimension  of  the  microstructural  features  in  which  cracks  initiate,  which  for  the  Ti-6246  alloy 
are  primary  alpha  particles  approximately  5  pm  in  diameter.  Under  the  latter  definition  of 
cracking,  our  consistent  results  on  repeated  tests  indicated  that  the  minimum  of  Nt  approached  0 
cycles,  which  meant  that  the  minimum  lifetime  was  spent  almost  completely  in  the  growth  of  a 
crack  that  begins  on  the  microstructural  scale.  Although  the  probability  of  occurrence  of  N  ->() 
may  be  reduced  by  factors  such  as  the  imposition  of  surface  residual  stresses  through  shot 
peening,  in  the  stress  range  of  rotor  design,  it  appears  difficult  to  preclude  the  rare  occurrence  of 
secondary  damage,  which  also  leads  to  Nt  ->0. 

Our  research  on  Ti-6246  and  other  high  performance  materials  suggests  that  the 
presumption  of  N,  ->0  for  rotor  design  stresses  may  have  a  level  of  general  applicability.  If  this 
is  correct,  it  may  be  attractive  to  reinterpret  conventional  approaches  and  strategies  for  life-cycle 
design  and  management  of  turbine-engine  rotor  materials  and  components. 
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Figure  1.  Schematic  of  design  requirements  and  expected  material/component  response 
specified  by  the  Department  of  Defense  Propulsion  System  Integrity  Program  (PSIP). 
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Figure  2.  Expanding  on  the  concept  of  “Expected  Material  Properties”  highlighted  at  the  upper 
left  of  Fig.  1,  it  is  commonly  assumed  that  fatigue  variability  is  well  described  by  a  single 
probability  distribution  (light  blue)  as  show  above.  Traditional  methods  of  design  and  life 
management  derive  material  life  limits  by  extrapolation  to  the  lower  tail  of  the  single 
distribution.  However,  it  now  appears  that  such  fatigue  may  actually  arise  from  at  least  two 
separate  underlying  probability  distributions.  In  this  sense,  the  identification  of  two  distributions 
appears  to  offer  substantial  opportunity  for  improved  understanding  of  the  fundamental  material 
behavior  and  a  reduction  in  predictive  uncertainty. 
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20pm 

Figure  3.  Nominal  duplex  microstructure  of  Ti-6Al-2Sn-4Zr-6Mo  alloy  processed  to  in  the  a+p 
phase  field,  containing  approximately  20%  primary  a  grains  and  fine  a-phase  platelets  in  a 


10  pm 

Figure  4.  Typical  fatigue  crack  origin  (lower  edge)  showing  a  single,  surface-connected,  alpha 
grain  facet  where  crack  nucleation  occurred. 
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Figure  5.  Ti-6246  stress  vs.  life  (S-N)  with  plot  showing  individual  specimen  lifetimes  and  a 
regression  line  representing  mean  behavior.  Tests  above  750  MPa  performed  at  20  Hz,  and  tests 
below  750  MPa  performed  at20  KHz. 


Cycles  to  Failure 

Figure  6.  Lognormal  probability  plot  of  fatigue  lifetime  data  for  individual  stress  levels  shown 
in  Fig.  5.  From  left  to  right:  amax  =1040,  925,  900,  860,  820,  700,  650,  600,  &  550  MPa. 
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Figure  7.  Lognormal  probability  plot  showing  820  MPa  and  860  MPa  data  from  Fig.  5, 
highlighting  three  regions  of  the  bimodal  trend  (Type  I,  transition,  Type  II). 
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Figure  8.  Lognormal  probability  plot  of  Type  I  and  Type  II  fatigue  lifetime  data  alone  for 
individual  stress  levels  shown  in  Fig.  5.  From  left  to  right:  amax  =1040,  925,  900,  860,  820,  860, 
820,  700,  650,  600,  &  550  MPa.  Note  that  820  and  860  MPa  data  each  appear  twice,  since  these 
stress  levels  exhibited  both  Type  I  and  Type  II  behavior.  As  shown  in  Fig.  6,  data  from  700  and 
900  MPa  tests  also  exhibited  both  types  of  behavior,  but  to  a  lesser  extent  than  at  820  and  860 
MPa. 
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Figure  9.  Lognormal  probability  plot  of  lifetime  data  for  tests  with  amax  =  820  MPa,  showing  a 
regression  extrapolation  to  the  B0.1  (1/1000)  probability  of  failure,  along  with  95%  confidence 
bands  for  the  fit. 


Figure  10.  Lognormal  regression  of  the  Type  I,  shortest-life,  points  from  alone,  showing  that  the 
slope  of  the  regression  is  much  steeper  and  the  B0.1  extrapolated  life  is  much  greater  than  shown 
in  from  the  full  dataset  in  Fig.  8. 
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Figure  11.  Histograms  of  primary  a  particle  sizes  (ap)  sizes  (blue)  as  measured  from 
metallurgical  section  and  sizes  of  crack  initiation  regions  (red)  measured  as  projected  normal  to 
the  fracture  surface. 


Figure  12.  Plot  of  da/dN  vs.  AK  for  several  naturally  initiated  small  surface  cracks  measured  via 
acetate  replication  of  the  specimen  surface,  and  regression  trend  from  similar  data  on  large 
cracks  in  conventional  compact  type  C(T)  specimens. 
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Figure  13.  Probability  plot  of  Monte  Carlo  model  simulation  of  Type  I  fatigue  lifetimes  based 
only  on  microstructure  (Fig.  12)  and  small  +  large  crack  (Fig.  13)  results  plotted  vs.  the  short  life 
data  from  the  from  the  amax  =820  MPa  experiment  (Fig.  9). 
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Figure  14.  S-N  plot  showing  (i)  individual  regression  lines  (blue  and  dashed)  to  Type  I  and  Type 
II  data,  (ii)  schematic  bimodal  probability  distribution  transitioning  from  Type  I  to  Type  II 
dominance  with  decreasing  stress,  and  (iii)  Monte  Carlo  predictions  of  50%,  10%,  and  0.1% 
probabilities  of  failure  for  the  minimum  lifetimes. 
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Figure  15.  Probability  of  occurrence  of  Type  I  (short-life)  fatigue  based  on  bimodal  regression 
analysis  of  fatigue  lifetime  data  from  Fig  14. 
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Figure  16.  (a)  1000  Monte  Carlo  model  predictions  of  crack  depth  vs.  life  (a-N)  for  a  maximum 
stress  of  820  MPa,  and  (b)  and  the  same  predictions  on  an  expanded  scale  that  shows  the  area  of 
primary  interest  for  design  and  life  management  of  rotor  components  designed  to  a  total  lifetime 
of  8000  cycles  (TACs) 
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